Hexagons become second if symmetry is broken 
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Pattern formation on the free surface of a magnetic fluid subjected to a magnetic field is inves- 
tigated experimentally. By tilting the magnetic field the symmetry can be broken in a controllable 
manner. When increasing the amplitude of the tilted field, the flat surface gives way to liquid 
ridges. A further increase results in a hysteretic transition to a pattern of stretched hexagons. The 
instabilities are detected by means of a linear array of magnetic hall sensors and compared with 
theoretical predictions. 

PACS numbers: 05.45-a, 47.54+r, 47.52 



Pattern formation in isotropic systems is more compli- 
cated than in anisotropic ones: One of the hallmarks 
of isotropic systems is the possibility to bifurcate to 
hexagons from an unstructured ground state, which is 
due to the existence and interaction of 3 degenerate wave 
numbers 1]. This situation is structurally unstable, how- 
ever: The smallest distortion of this symmetry acts as 
a singular perturbation and will lead to a qualitatively 
different instability, namely a primary bifurcation to a 
stripe like pattern. A specific example has recently been 
calculated in detail |2| for a magnetic fluid |3(. In the 
ideal isotropic system, hexagons will occur under the in- 
fluence of a magnetic field which is perfectly normal with 
respect to the fluid surface 4] . The slightest change of the 
orientation of the magnetic field is predicted to change 
this subcritical transition: Ridges appear supercritically 
via the primary bifurcation. Their interaction with waves 
along the less- favoured direction gives rise to " stretched" 
hexagons via a secondary bifurcation. 
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FIG. 1: Surface patterns of magnetic fluid in a magnetic field 
tilted by the angle ip — 23° to the vertical. Surface reflections 
of the liquid ridges for the verical induction B =20 mT (a) 
and of the tilted crests for B =32 mT (b). The side view of 
the patterns is presented in (c) and (d), respectively. 



A first observation of liquid ridges was reported in Ref. 
|3- In this paper we present a quantitative characteri- 
zation of the primary bifurcation to liquid ridges and a 
secondary bifurcation to a pattern of stretched hexagons, 
as shown in FigQJ via use of a magnetic measurement 
technique. Specifically, we measure the threshold induc- 
tion B p and B s for the primary and secondary instabil- 
ity for various angles of tilt (p. The measurements of 
B p (cp) agree with the theoretical prediction if the non- 
linear magnetization curve of the magnetic fluid used in 
the experiment is taken into account. 
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FIG. 2: Scheme of the experimental setup. 



Our experimental setup is shown in Fig.EI A cylin- 
drical vessel with an edge machined from Teflon® with 
a diameter of 12 cm and a depth of 2 mm is brimful 
filled with fluid and is situated in the center of a pair 
of Helmholtz coils. For details see Ref. 6]. The axis of 
the coils can be tilted against the vertical by an angle 
(p = [0°,90 ]. The experiments are performed with the 
magnetic fluid EMG 909 Lot F061998B (Ferrotec Corp.), 
with fjL r = 2.11. A charge-coupled-device (CCD) camera 
is recording the patterns from above. In order to measure 
the amplitude of the steep crests a linear array of 32 hall 
sensors was mounted 1.78 ± 0.1 mm below the bottom of 
the dish. By this technique the local increase of the mag- 
netic induction below a liquid crest is utilized to measure 
its amplitude. The sensors communicate via 32 amplifiers 
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and a bus with the PC. Details of this method are pre- 
sented elsewere |g • For calibration purpose a commercial 
Hall probe (Group3-LPT-231) in combination with the 
digital teslameter (DTM 141) was used. 

The magnetic field is tilted towards the x-axis. Increas- 
ing the magnetic induction, we observe a transition from 
the flat layer of magnetic fluid to the pattern of liquid 
ridges, displayed in FigCJi,c. The wave vector of the pat- 
tern is oriented along the y-axis and thus perpendicular 
to the horizontal field component. The vertical compo- 
nent of the local magnetic induction was measured by 
means of the sensor array oriented parallel to the wave 
vector. In order to reduce the spatial inhomogeneities 
of the magnetization caused by the finite container size, 
the spatial variation measured at a subcritical induction 
of 20.5 mT is substracted. The ensuing local induction 
B(x) is presented in Fig.Olfor different values of the ap- 
plied magnetic field, measured at the tilt angle cp = 32°. 
The open circles mark the data, the solid line gives the 
least square fit to 



B(x) = A • cos(kx -i/j) + B. 



(1) 



Here A denotes the modulation amplitude, k the absolute 
value of the wave vector, xp the phase, and B the mean 
value of the induction. 
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FIG. 3: Profiles of B(x), for different values of the applied 
magnetic field. The open circles mark the data, the solid lines 
the fits by Eq. Q. 



The square of the modulation amplitude A is plotted in 
Fig.^ versus the control parameter B. The monotonous 
increase after a threshold B p is characteristic for a super- 
critical bifurcation. It can be described by the solution 
of the stationary amplitude equation Q 



= e p A - gA 3 + b. 



(2) 



In accordance with the symmetry of the problem 
e p = (B 2 — B 2 )/B 2 was selected to be the dimensionless 
bifurcation parameter, g is the cubic coefficient a scaling- 
, and b a imperfection parameter. The solid line in Fig^I 
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FIG. 4: Square of the modulation amplitude of the liquid 
ridges versus the mean magnetic induction B of the array 
detector. For details see text. 



gives the fit of the experimental data by the solution of 
Eq.©. We obtain B p = 21.17mT, g = 21.16mT- 2 , and 
the slight imperfection b = 4.3 • 10 _5 mT. The dotted line 
displays the solution without imperfection (b = OmT). 

Increasing the controll parameter further initiates a 
secondary instability to the stretched hexagonal pattern 
as displayed in Fig|H),d. The blow-up in (d) indicates 
that the crests riding on top of the ridges are asymmet- 
ric with respect to the wave vector of the ridges. Thus 
the pattern (of stretched hexagons) lacks any non-trivial 
rotational symmetry. 

For a quantitative analysis of the secondary instability 
a series of 400 measurements of the local induction B(x) 
has been performed for cp = 23°. For clarity Fig. [SJi (b) is 
presenting only every 20th line for a quasistatic increase 
(decrease) of the control parameter 5, respectively. In 
order to detect both, the ridges and the crests, the sensor 
line is now oriented with an angle uo — (75.8 ± 0.05)° to 
the y-axis. In this way it is covering 2 | ridges, which 
can be recognized in the lower part of the plots. For B « 
22 mT the transition to the stretched hexagons occurs. 

A mathematical characterization of the stretched 
hexagons can be obtained as follows. In a stretched- 
hexagonal pattern the wave vectors fulfill the side con- 
dition — k\ = &2 + ks and k = = \ks\. With the 

k/k and b 



abbreviations k = |fei|, n 



= V4n 2 - 1 the 
wave vectors read k\ = fc(0, 1,0), k<i = — 1(— 6, 1,0), and 
ks = — 1(6, 1,0) which coincide for n = 1 with the vec- 
tors for a regular hexagonal pattern. The amplitude of 
the ridges A^ and of the stretched-hexagonal pattern Aq 
can be combined to the amplitude of the over-all pattern 



A(x) Aq cos kix + 22 
3 i=i 



COS kj2 



(3) 
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FIG. 5: Local magnetic induction for increasing (a) and 
decreasing (b) control parameter B and ip = 23°. The solid 
lines give the fit by Eq. (0. 



A cut through this pattern is given by 

xm = x o + t • e (4) 

where t denotes the distance from the starting point xq = 
(xo,2/o 5 0) of the sensor line and e= (sin a;, cos a;, 0) its 
unity vector of orientation. Plugging Eq.(@J) in Eq.© 
yields 

A(x M ) = A(t) = A*R(t) + A*H(t). (5) 
for the amplitude along the sensor line. Here 

R(t) = cos (k(yo + t coscj)) (6) 
gives the contribution of the ridges, and 

H(t) = - [cos(k(y + tcosu)) + cos(fe(^* + 

+ cos(fc($ + (7) 

the contribution of the hexagons. Here = 
h(bx +yo) , ^ = ^(bsmco + cosco) , = 
2 (foco — ) and \£* = \ (b sin a; — cos a;) are abbrevi- 
ations. For small amplitudes Eq.Q is sufficient, but for 
higher amplitudes it is important to take into account 
the higher harmonics k m = m • k with rn = 1,2, ... The 
surface is than given by 

Mr M H 

A(t) = £ A^R m (t) + £ A^H m (t). (8) 

771=1 171=1 



This model is fitted to the data by four nonlinear pa- 
rameters, which are the wave number k of the ridges, the 
starting point x~o = (xo,yo) of the sensor line, and the 
stretching factor n of the hexagonal pattern. The am- 
plitudes Aq and Aq are linear parameter of the basic 
functions R(t) and H(t). The solid lines in Fig|^|give the 
best fit by Eq.® taking into account the basic mode of 
the ridges (Mr = 1) and the first two of the hexagons 
(M H = 2). 

From this fit the amplitude A$ of the hexagons can 
be extracted. It is plotted in Fig|£|vs. the control pa- 
rameter B. The open squares (circles) mark the data for 
an increase (decrease) of 5, respectively. The hysteresis 
is characteristic for a subcritical bifurcation, which has 
been predicted for the transition from ridges to stretched 
hexagons Q. 
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FIG. 6: Amplitude A of the stretched hexagonal pattern 
versus the applied magnetic induction. The open squares 
(circles) denote the measured amplitude under increase (de- 
crease) of the induction, respectively. For clarity only ev- 
ery 10th point is shown. The solid (dotted) line give the 
fit by Eq.© with b = 1.6 • 10" 4 mT (b = OmT), respec- 
tively. For the other parameters we obtained B s = 22.594mT, 
7l 7.6mT"\ 72 0.9mT"\ and g 116.31mT" 2 . 

Next we describe the amplitude Ah = Aq of 
the hexagons after the secondary bifurcation at 
es = (B 2 — B$)/Bp. In the spirit of a weakly nonlinear 
analysis slightly above es we use the amplitude equation 

= eA H + 7i (1 + j2es)A 2 H - gA% + b s . (9) 

In this experimental paper the coefficients in Eq.® have 
been obtained by a fit to the measurements in order to 
circumvent their tedious calculation from the basic equa- 
tions. To avoid the ambiguity of Au{B) in the hysteretic 
regime, B(Ah) was fitted to the data according to Ref. 
Q. The result of the fit is presented in Fig|£|by a solid 
line, while the dashed line gives the solution for the same 
parameters, however with bs = 0. 
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For decreasing B the system follows the solid line very 
well down to the saddle- node. For increasing B the agree- 
ment is less convincing in the bistable regime. Here the 
impact of the edges [lOj seems to penetrate the interior 
of the dish much stronger. As a consequence the analysis 
by Eq.Q is not sufficient in this regime - see also Fig|^- 
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FIG. 7: Critical inductions versus the inclination of the mag- 
netic field. The data marked by open circles have been es- 
timated by fitting the evolution of the ridge amplitude (e.g. 
Fig.^J by Eq. @. The other values, marked by full symbols, 
have been determined by visual inspection of the liquid layer. 
The solid line gives the fit obtained via Eq. (|10f> . 



To conclude, for the tilted field instability we have mea- 
sured the forward bifurcation to liquid ridges. The an- 
gular dependence observed in the experiment, is quanti- 
tatively described by taking into account the nonlinear 
magnetization. In addition we measured the backward 
bifurcation to hexagons, which has been predicted by an 
energy variational method . A full quantitative agree- 
ment with these predictions can not be expected, be- 
cause the theory is restricted to permeabilities \i r < 1.4, 
while we had to use \i r = 2.11 to avoid huge fields. The 
essence of the experimental observation, namely a struc- 
tural change of the primary instability, seems to be well 
described by this theoretical ansatz: For broken symme- 
try ridges always precede hexagons. They are increas- 
ingly difficult to resolve, however, if the angle of tilt di- 
minishes. A similar scenario can be expected e.g., for 
non-boussinesq inclined layer convection |l4|, for mag- 
netohydrodynamic as well as electro-convection in tilted 
magnetic fields [l5(, for lucent hexagons under influence 
of an asvmmetric Fourier filter and for Turing pat- 
terns |l6| in stressed gel. 
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Next we investigate the angular dependence of the crit- 
ical induction for the first and secondary bifurcation. In 
Fig.0the measured data for the transition to ridges are 
marked by circles, whereas the transition to stretched 
hexagons is denoted by triangles and the reverse transi- 
tion to ridges by squares. 

The solid line gives B p calculated for the instability 
of the fiat surface. It follows from the dispersion rela- 
tion of the surface waves in ^/-direction with uo 2 = and 
k y = k c = y/pg/cr, using the relation between magnetic 
induction and magnetization, 

U\k x =^ky) = g\ ky \ + -\ky\ 3 -^k 2 y M 2 Z (10) 



M 2 - M 2 ^* + M J M 2_ M 2^x' 

z x+i x+i V z x +1 

The dispersion relation Ea.(|lUj) takes into account the 
nonlinearity of the magnetization curve M(H) of the 
fluid and can be deduced from Eq.(36) in 11] whereby 
M = Mz& z + Mx& x denotes the magnetization of the 
fluid for the undisturbed surface. The susceptibilities 
x(H) = M ^ and x\H) = dM Q ^ were determined 
from the experimental magnetization curve assuming a 
logarithmic normal distribution for the size of the mag- 
netic particles in the fluid 05 0. In contrast, in Ref. Q 
a constant x nas been used, which results in a B p not 
depending on tp. 
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